J/^Ij suppression in ultrarelativistic nuclear collisions 



Bin Zhang*^, C. M. Ko"", Bao-An Li^, Ziwei Lin^, and Ben-Hao Sa'' * 
^Cyclotron Institute and Physics Department, 
Texas A&M University, College Station, TX 77843, USA 
^Department of Chemistry and Physics 
Arkansas State University, P.O. Box 4^9 
State University, AR724 67-0419, USA 



Abstract 



Using a multiphase transport model, we study the relative importance of 
J/ip suppression mechanisms due to plasma screening, gluon scattering, and 
hadron absorption in heavy ion collisions at the Relativistic Heavy Ion Col- 
lider. We find that for collisions between heavy nuclei such as Au+Au, both 
plasma screening and gluon scattering are important. As a result, the effect 
due to absorption by hadrons becomes relatively minor. The final J/ij: survival 
probability in these collisions is only a few percent. In the case of collisions 
between light nuclei such as S+S, the effect of plasma screening is, however, 
negligible in spite of the initial high parton density. The final J /ip survival 
probability thus remains appreciable after comparable absorption effects due 
to gluons and hadrons. 
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I. INTRODUCTION 



According to the fundamental theory of the strong interactions, Quantum Chromo- 
Dynamics (QCD), normal nuclear matter is expected to undergo a phase transition to de- 
confined quarks and gluons when its density and/or temperature are high To produce 
such a quark-gluon plasma (QGP) in the laboratory, experiments involving collisions of nu- 
clei at relativistic energies have been carried out at the CERN SPS 0. Possible evidences 
for the production of this new phase of matter have recently been announced 0. Further 
experiments at the Relativistic Heavy Ion Collider (RHIC), which allows collisions at much 
high energies than those available previously, are expected to provide a better opportunity 
to create the quark-gluon plasma and to study its properties. 

Since quarks and gluons cannot be directly detected in experiments, many indirect ob- 
servables have been proposed as possible signatures for the QGP. These include: enhanced 
production of strange hadrons as a result of the short strangeness equilibration time in QGP 
1^, suppression of J /ip production due to color screening in QGP quenching of high pt 
jets due to passage through the QGP and enhancement of low mass dileptons as a 

result of medium modifications of hadron properties in high density matter P,p!0[]. Although 
all these signals have been observed in heavy ion collisions at CERN SPS, alternative expla- 
nations without invoking the formation of the quark-gluon plasma have also been proposed. 
As the QGP is expected to be produced at RHIC, it is of interest to study these signatures 
in heavy ion collisions at such high energies. Since the quark-gluon plasma has a finite size, 
exists for a finite time, and may not be in equilibrium, it is important to use a dynamical 
model to take into account these effects. Using the parton cascade model Satz and 



Srivastava [12| have recently studied the time evolution of the density profile of the parton 
system in these collisions and demonstrated the possibility of using a dynamical model to 
study the onset of deconfinement and its effect on quarkonium suppression. In this paper, 
we shall investigate explicitly J / ip suppression in central nucleus-nucleus collisions at RHIC 
energies using a multiphase transport model (AMPT) |]13[ . Since the AMPT model includes 
both initial partonic and final hadronic interactions as well as the transition between these 
two phases of matter, it allows us to study the relative importance of the partonic and 
hadronic effects on J/ip suppression. 

This paper is organized as follows. In Section we briefly describe the AMPT model 
for ultrarelativistic heavy ion collisions and its extension to include J/ip production. The 
mechanisms for J /if) suppression are given in Section IH. Results for Au+Au and S-l-S 
collisions at RHIC are given in Section |IV|. First, we show the time evolution of the particle 
number and energy densities in both the partonic and hadronic matters. This is then followed 
by the time evolution of various J/ip suppression effects and the J /if) survival probability. 
Finally, a summary is given in Section ^ 



II. MULTIPHASE TRANSPORT MODEL 

In the AMPT model, the initial conditions are obtained from the minijets generated by 
the Heavy Ion Jet Interaction Generator (HIJING) [l^ by using a Woods-Saxon radial shape 
for the colliding nuclei and including the nuclear shadowing effect on partons via the gluon 



recombination mechanism by Mueller-Qiu |]T5|. After the passage of the colliding nuclei. 
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the Gyulassy-Wang model is then used to generate the initial space-time information of 
partons. The subsequent time evolution of the parton phase-space distribution is modeled by 
the Zhang's Parton Cascade (ZPC) |T^,[T^, which at present includes only the gluon elastic 
scattering. After partons stop interacting, the Lund string fragmentation model [p!8| , p!9| is 
used to convert them into hadrons. The dynamics of resulting hadronic matter is described 
by a relativistic transport model (ART) [^]. Details of the AMPT model can be found in 

To include J/ip in the AMPT model, we use the perturbative approach [|1^ as its pro- 
duction probability is small in heavy ion collisions. Specifically, the reaction gg —>■ J/ijjg 
is selected from the PYTHIA program whenever there is an inelastic scattering between 
projectile and target nucleons. Instead of allowing the produced J/ip to undergo multiple 
interactions, fragmentation, and decay in the PYTHIA, its momentum is stored in a file 
and read into HIJING. The transverse position of the J/ip is determined by propagating it 
to the time of formation from the average transverse position of the colliding projectile and 
target nucleons. Since the J/ip production probability is increased from a{NN J/ipX) to 
ainei, each J/ip is given a probability of a{NN J /ipX) / ainei in order to have the correct 
number in an event, given by the total probability of all J/ip^s. The J/ip formation time is 
taken to be 0.5 fm/c, which is suggested by the virtuality argument of Kharzeev and Satz 
pT| and is also consistent with the uncertainty principle used for estimating the lifetime of 



an expanding cc [E2],E3l. Before this time, the cc pair is considered a precursor J/tp. 



III. J/ij ABSORPTION 

After a pair of cc is produced in the initial collisions, whether it can materialize as a, J/ifj 
depends on its interactions in the initial partonic and final hadronic matters. In the partonic 
matter, the cc pair, which are initially close in phase space and would normally form a J/^ 
after the formation time of 0.5 fm/c, will move beyond the confinement distance of J/ip and 
become unbound if the plasma screening remains strong [Q. They will later combine with 
the more abundant light quarks to form instead the charm hadrons when the parton density 
is low. The critical density for the plasma screening to be effective can be estimated using 
ric = (2ko^1) / {?>n'^as) 101 5 where /Xc is the critical Debye screening mass, as is the QCD 
coupling constant, and fco is the slope parameter of the transverse momentum distribution. 
Using lie = 0.7 GeV, as = 0.47, and = 0.6 GeV, one obtains a critical density of 5 
fm~^. Following the method of Ref. |[T^, we first determine from the AMPT model the 



time evolution of the radius of the volume in which the parton density is above the critical 
density. A J/ip is thus not formed from the cc pair if its radial position after the formation 
time of 0.5 fm/c is within the critical radius at that time. 

Besides dissociation due to plasma screening, a. J/ip may also be dissociated by collisions 
with gluons in the partonic matter [^. The cross section for gJ /ip cc has been estimated 



in Ref. p5[ to be 3 mb. Similarly, a precursor J/ip can also be dissociated by gluon scattering. 
The cross section is expected to be smaller as its size is smaller than a physical J In the 
present study, we ignore the difference between the scattering of a gluon with precursor and 
physical J/ip., and use the same cross section. Since this cross section is not well-known, we 
have also used a value of 1 mb in the following study. 

In the hadronic matter following the parton stage, a, J/ip can be further destroyed by 
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FIG. 1. Time evolution of particle number and energy densities in the central cell around 
r] = and r = for central (b=0) Au+Au collisions at RHIC. 



collisions with hadrons. The J/ip absorption cross section by baryons is taken to be 6 mb 
26| while that by meson is taken to be 3 mb ||2^ above their respective thresholds. Although 



earlier studies based on the perturbative QCD p8[ and simple hadronic model |29 



much smaller J/ijj absorption cross sections by hadrons, the above values are consistent 



with recent studies using the quark-interchange model and the more complete hadronic 
model [^l],^. Using these cross sections in the transport model, it has been found that the 
observed suppression of J/ip production in heavy ion collisions at the SPS energies can be 
reasonably described except for the most central collisions ||33|-|35| . As the formation time of 



J/ip is comparable to the passing time (0.1 fm/c) of two colliding nuclei at RHIC energies, 
J/ip absorption by the projectile and target nucleons is not expected to be important ||3^ 



and will be ignored in the present study. We note that as the J/ ijj is treated perturbatively 
in the simulation, the momenta of other hadrons are not affected by their interactions with 



IV. RESULTS 

A. time evolution of particle number and energy densities 

In Fig. |l|, we show the time evolution of the particle number and energy densities of 
partons and hadrons in central Au+Au collisions at a/s = 200 GeV. They are the densities 
in the central cell, which has a transverse radius of 1 fm for partons and 2 fm for hadrons. 
The longitudinal dimension of the central cell is taken to be ±5% of the time t, which is 
equivalent to taking the central space-time rapidity cell. The results are not significantly 
changed if ±10% of the time is used for the longitudinal dimension. We see that the initial 
parton density is about 30 fm~^ and is much higher than the J/ ip dissociation critical density 
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FIG. 2. Same as Fig. |T| for central (b=0) S+S collisions at RHIC. 



(about 5 fm~^) given above. The time evolution of these densities is seen to deviate from that 
based on the ideal Bjorken boost invariant scenario [38|, which would lead to a linear curve 
on the log- log plot. This difference is mainly due to the more realistic treatment of initial 
collisions by using the the Gyulassy-Wang model for the formation time and the presence 
of radial flow as well as a gradual freeze-out |^ at the later hadronic stage. From the ratio 
of the parton energy density to its number density, one sees that the average energy of a 
parton is more than 1 GeV. The parton stage lasts about 2 — 3 fm/c, while the hadron stage 
starts gradually at around 3 — 4 fm/c and lasts until about 10 fm/c. The average energy 
of a hadron is less than 1 GeV, since the central rapidity is dominated by mesons instead 
of baryons in ultrarelativistic nuclear collisions. In Fig. |], we show the results for central 
collisions of S+S at ^/s = 200 AGeV. Due to the smaller size of the system comparing to 
that of Au+Au collisions, the plasma lifetime is shorter, i.e., about 1 — 2 fm/c, and the 
hadron stage sets in at about 2 — 3 fm/c. We note that even in the smaller parton system 
produced in S+S collisions, the initial density is much higher than the critical density for 
J/ip dissociation. 

In Fig. 1^ and Fig. ^, we show, respectively, the parton density in central Au+Au and 
S+S collisions at different times. As minijet gluons are produced from initial hard collisions, 
their densities at different radii reflect the number of initial binary nucleon-nucleon collisions. 
As seen from the figure, they are different from the initial nuclear density distribution given 
by a Woods-Saxon form. As expected, both the size and lifetime of the partonic matter 
produced in S+S collisions are smaller than those in Au+Au collisions. From the time 
evolution of the parton density, one can determine the time evolution of the critical radius 
for J/ ip dissociation, and this is shown in Fig. ^ by full dots for Au+Au collisions and open 
dots for S+S collisions. The solid lines are polynomial fits to the above results using the 
form, rc(t) = flo + ^it + a2t^ + a^t^ for tmin < t < tmax, with rdt) = rc{tmin) 

for t <^ imin 

and rc(t) = for t > tmax- For the partonic matter produced in Au+Au collisions, we have 



0.15 fm/c, U 



1.2 fm/c, Oo = 6.39, ai = —4.06, 02 = 4.79, and = —4. 
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FIG. 3. Parton density in central (b=0) Au+Au collisions at RHIC at different times. The 
numbers are the proper times in fm/c for the density measurements. 
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FIG. 5. Time evolution of the critical radius for plasma dissociation of J/ij: in central Au+Au 
and S+S collisions at ^/s = 200 AGeV. 

For S+S collisions, they are tmin = 0.12 fm/c, tmax = 0.5 fm/c, ao = 3.278, ai = —11.949, 
02 = 31.514, and 03 = —41.428. In Au+Au collisions, the critical radius for J/tp dissociation 
extracted from our model is about 6 fm at beginning of the parton cascade and vanishes 
after about 1.2 fm/c. The critical radius is reduced to about 2 fm initially and vanishes 
after only about 0.5 fm in S+S collisions. Since the duration of the partonic matter that is 
above the critical density for J/ip dissociation is shorter than the J/ip formation time, J/i/) 
suppression due to plasma screening is therefore unimportant in central S+S collisions. 

B. time evolution of J ftp survival probability 

Including the above three mechanisms for J/ip suppression in the AMPT model, we 
have evaluated the J/ ip survival probability in ultrarelativistic heavy ion collisions at RHIC 
energies. In Fig. ^ we show for central Au+Au collisions the time evolution of the J/ijj 
formation probability Pf(t), the J/ip absorption probability P^it) by gluon scattering, the 
J/ip dissociation probability Pdit) by plasma screening, the J/ip absorption probability Pcit) 
by hadrons, and the J/ip survival probability Psit) = 1 — -P|(t) — Pd(t) — Pc(f)- The results 
are taken for the J/ip produced in the central rapidity interval < 1. It is seen that 

the J/ip formation probability increases quickly with time and about 90% of the J/ip are 
formed by 0.5 fm/c. Absorption by gluons starts very early in the process and ends at about 
1 fm/c after dissociation due to plasma screening begins. The latter also ends at about 1 
fm/c. Both dissociation due to plasma screening and absorption by gluon collisions give 
comparable contributions, accounting for the suppression of about 90% of the J/ip, while 
absorption by hadrons contributes only about a few percent. The final J/ip survivability is 
about 6%. 

The results for the case without plasma screening is shown in Fig. ^. Comparing to Fig. 
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FIG. 6. Time evolution of the J /ip formation probability Pf{t), the J/ip absorption probability 
P^{t) by gluon scattering, the J/^p dissociation probability Pdit) by plasma screening, the J/^p 
absorption probability P^{t) by hadrons, and the J ftp survival probability Ps{t). For the survival 
probability, the solid curve includes absorption by both gluons and hadrons, while the long-dashed 
curve includes only absorption by gluons. 

^, we see that the J/ip absorption probability by gluons is similar for time up to 0.5 fm/c, 
but it lasts longer until 2 — 3 fm/c. This indicates that there is a competition between J/ijj 
suppression due to gluon scattering and plasma screening, i.e., some J/V^'s that are destroyed 
by collisions with gluons would have been dissociated by plasma screening if they are not 
allowed to scatter with gluons. Fig. |] shows the results without J/ ^p absorption by gluon 
scattering. In this case, the plasma dissociation probability saturates quickly. Compared 
to that shown in Fig. ^, the J/ip suppression probability is seen to start late but saturate 
early. This difference in the survival probability may be seen by studying the azimuthal 
distribution of final survival J/ ip in mid-central collisions. It is interesting to note that J/ i/j 
suppression during the parton stage is similar whether when both plasma screening and 
gluon scattering are present or when only one of them is present. The effect of reducing the 
J/ ip scattering cross section by gluon is shown in Fig. ^. Compared with Fig. we see that 
the decrease in the contribution from gluon scattering is partly compensated by the plasma 
screening. 

In Table 0, we summarize the final absorption probabilities due to gluons (P^), hadrons 
(P^^), plasma dissociation (Pd), and the survival probability (Pg) in central Au+Au collisions. 
It is seen that the J/ip survival probability after the parton stage, Sg = 1 — Pc — Pd, is about 
9% and 17% for Cgj/^ = 3 mb and 1 mb, respectively. The parton stage thus has a large 
effect on J/ip suppression in heavy ion collisions at RHIC. In both cases, the J/ip suppression 
factor in the hadron stage is Sh = Phi Sg ^ 35%. Depending on the different assumptions on 
the mechanisms for J /ip suppression, the final survival probability can range from about 6% 
to about 26%. This indicates that finite size effects may be important in J /ip suppression at 
RHIC energies. We have also studied J/ip suppression in S+S collisions at RHIC energies, 



Au(1 00AGeV)+Au(1 OGAGeV) 
(b=0,o„j,^=3mb,|yj,^|<1) 
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FIG. 7. Same as H without plasma dissociation. 
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FIG. 8. Same as |^ without gluon destruction. 
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FIG. 9. Same as ^ with reduced J absorption cross section by gluons. 

TABLE I. Final J /ip absorption probability P§ by gluon scattering, J /ip dissociation probabil- 
ity Pd by plasma screening, J /ip absorption probability by hadrons, and J/tp survival probability 
Pg for Au+Au (6 = 0, y/s = 200 AGeV). {gJ/ip)a indicates cr^j/^ = 3 mb while {gJ/'ip)b indicates 



o-gj/t/. = 1 mb. 

Pj Pd Pj' Ps 

screening+(c/J/V')a+/iJ/V' 43.9% 47.1% 2.9% 6.1% 

screening+(5J/V')fe+/iJ/'0 20.6% 62.4% 6.5% 10.5% 

screening+hJ/V' 75.4% 10.2% 14.4% 

{gJ/'ip)a+hJ/iJ 78.3% 10.6% 11.0% 

lgJ/i!)b+hJ/i) 42.9% 31.1% 26.0% 

hJ/^ 56.7% 43.3% 



and the results are summarized in Table |T[ As expected from the discussion of Fig. ^ there 
is no contribution from plasma screening in collisions of such light system, although the 
initial density is above the critical density. As a result, the J /if) has an appreciable survival 
probability (about 50%) after absorption by gluons and hadrons. Since J /if) absorption by 
gluon scattering still contributes appreciably to the final J/ip suppression, it may provide 
an opportunity for studying this mechanism. We note that the J /if) suppression factor in 
the hadron stage is about Sh ~ 16%, which is smaller than in Au+Au collisions. 

V. SUMMARY 

In summary, we have studied J/ip suppression in Au+Au and S+S collisions at RHIC 
energies in a multiphase transport model. It is found that finite size effects may be im- 
portant in J/ifj suppression at RHIC energies even in central Au+Au collisions. Both J/ip 
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TABLE II. Same as Table | for S+S {b = 0,^ = 200 AGeV). 
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suppression due to plasma screening and gluon scatterings are important in Au+Au colli- 
sions, leading to a final J/ip survival probability of only about 6%. For S+S collisions, even 
though the initial energy density is above the critical density for J/ip dissociation, only J/ijj 
absorption by gluon and hadron scatterings contribute to its suppression as a result of the 
short lifetime of the plasma. J/ip suppression in heavy ion collisions between light nuclei 
thus provides a possible opportunity for studying the interactions between J/ and gluon. 

In our study, we have ignored both J/ ip formation from the recombination of cc pairs dur- 
ing hadronization and J/ ip production from charm mesons in the final hadronic matter. The 
number of J/ip formed from charm quarks during hadronization can be roughly estimated 
from considerations of the phase space and the spin and color factors. At the end of parton 
stage, the volume of the partonic matter in the central unit of rapidity is approximately 
given hj V = ttR'^t, where R and r are the transverse radius and longitudinal dimension, 
respectively. For Au+Au collisions, they are i? ~ 7 fm and r ~ 3 fm. The probability that 
an anticharm quark is found within a distance of the J/ ip radius from a charm quark is given 
by the ratio of the size of the J/ip to the size of the partonic matter. With a. J/ip radius 
^■j/i> ~ 0.3 fm, one obtains the probability 2.5 x 10~^. In order for this charm pair to form 
a J/ipi their spins must be parallel and their total color should be neutral. This reduces the 
probability by the factor (3/4) ■ (1/9) ~ 0.083. Requiring that the charm pair has a relative 
momentum distribution within po ~ l/'^J/V' ~ 0-67 GeV at T ~ 200 MeV further suppresss 
the probability by a factor of 0.29. Since there are about 1.5 charm pairs per unit rapidity 
per central Au+Au collision at RHIC, the number of J/ip formed from their recombination 
is about 1.5^ ■ 2.5 x 10~^ ■ 0.083 • 0.29 ~ 1.4 x 10~^, which is much smaller than the expected 
number of J/ip., i.e., 1.5/40 ■ 0.061 ~ 2.3 x 10^^. In the above, the factor 40 is the ratio of 
charm to J/ip production in initial hard collisions and the factor 0.061 is the final J/ip 
survival probability from Table |. As to J/ip production from charm mesons in hadronic 
matter, it has been shown to be insignificant in heavy ion collisions at RHIC, although it 



may not be negligible at LHC PT|,^ 



We have not included the effect due to ip' and Xc which are also produced in initial 
collisions and can decay to J/ip. Since these particles are less bound than J/ip, they are 
more likely to be dissociated and absorbed in both the quark-gluon plasma and the hadronic 
matter. As a result, inclusion of ip' and Xc is expected to lead to a smaller number of J/ip 
than predicted in the present study. Also, we have used the default gluon shadowing in 
HIJING without flavor and scale dependence, which is different from the one suggested by 
Eskola H3] that includes such dependence. If the latter is used, then we expect a reduction 



of the J/ip survival probability as the gluon density would be higher at RHIC energies [44 
These effects will be studied in the future. 
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